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ABSTRACT 

In this paper we use coherently integrated visibilities (see separate paper in these proceedings^) to measure the 
properties of binary stars. We use only the phase of the complex visibility and not the amplitude. The reason 
for this is that amplitudes suffer from the calibration effect (the same for coherent and incoherent averages) and 
thus effectively provide lower accuracy measurements. We demonstrate that the baseline phase alone can be used 
to measure the separation, orientation and brightness ratio of a binary star, as a function of wavelength. 



1. INTRODUCTION 

Binary stars are important for calibrating evolutionary 
stellar models. Because stellar models are sensitive to 
the parameters of the stars, it is important to obtain the 
highest possible accuracy of the stellar parameters. In 
this paper we will demonstrate how to extract bright- 
ness ratio and separation vector from measurements of 
the complex visibility phase only. This is significant be- 
cause the visibility phase does not suffer from the cal- 
ibration effects of visibility amplitudes, and therefore 
much higher precision can be obtained. Further, there 
are more complex visibility baseline phases than clo- 
sure phases such that using baseline phases instead of 
closure phases yields more information. Finally, com- 
plex visibility phases generally have better SNR than 
closure phases. 

2. THEORY 

The complex visibility of a binary star is 



V = cos (rj) + r cos (7) + i [sin (r"f) + r sin (7)] 
where 



7 = 



2tiB-s 
{r + l)X' 



B is the baseline vector, s is the separation vector, r 
is the brightness ratio, and A is the wavelength. The 
visibility phase is 



9 = tan 



with uncertainty 



sin (r7) -|- r sin (7) 
cos ir^) -\- r cos (7) 
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where N is the total number of photons counted and V 
is the coherently integrated visibility amplitude. 

In addition to the source phase, the measured phase 
also contains an instrumental phase component and a 
combination of atmospheric and vacuum phase terms. 
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The instrumental phase term can be measured by ob- 
serving a calibrator star (which has zero source phase) . 
The observed phase of a calibrator may contain both in- 
strumental and atmospheric phase terms, which is not 
a problem because the phase terms are additive. Fig- 
ure 1 shows the measured phase of a calibrator star. 
The atmospheric phase term takes the form 



2'K 

9atm = ^ [w + (n - 1) a] 
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Figure 1. The phase of a caUbrator star. This is the sum of the 
instrumental phase and some atmospheric phase. 




Figure 2. Example of what a binary star baseline might look 
like in a fringe-tracking interferometer. 



where z; is a vacuum path, a is the atmosphere path, n 
is the wavelength-dependent index of refraction, and 
is a phase offset. 

We can then measure the brightness ratio and sep- 
aration vector by fitting equation 2 to the basehne 
phases, with equation 1 inserted for the source phase, 
and equation 3 inserted for the atmospheric phase term. 
Note that the binary phase term in equation 1 is mono- 
tonic as it stands, but in a typical measurement we will 
see the phase be roughly centered around zero. This 
is because the fringe-tracker of the interferometer will 
add the appropriate vacuum delay to give the fringe 
close to zero phase. Figure 2 shows an example of what 
the phase of a binary might realistically look like in at 
high spectral resolution in a fringe-tracking instrument 
interferometer like the NPOI.^ 

3. DATA SETS 

We will present measurements of two stars, (fl Tauri 
and K Ursa Maj oris observed on the same night. On 
that night the NPOI was configured to observe three 
baselines, one of which was observed on two spectro- 
graphs. Table 1 lists the target observations. In addi- 
tion to target observations there were several calibra- 



Date 


Time (UT) 


Target 


Hour angle 


2004/1/30 


04:43 


^^Tau 


1.373 


2004/1/30 


06:09 


6*2 Tau 


2.810 


2004/1/30 


07:15 


6*2 Tau 


3.923 


2004/1/30 


06:37 


KUMa 


-1.284 


2004/1/30 


09:31 


KUMa 


1.612 


2004/1/30 


10;26 


KUMa 


2.499 


2004/1/30 


10:47 


KUMa 


2.871 


2004/1/30 


11:26 


KUMa 


3.541 



Table 1. Observations of d'^ Tauri. All observations lasted 30 s. 



Figure 3. UV coverage on d Tauri. 

tor observations to measure the instrumental phase for 
subtraction. 

4. ANALYSIS 

The coherently integrated visibilities are produced us- 
ing the methodology presented in the paper-'^ in these 
proceedings. The fit of the functional form of equation 2 
is not trivial. There are many false minimas such that 
it is important to either use a minimization method 
which finds the global minimum and ignores local min- 
ima, or to give the minimization procedure an initial 
guess which is within the capture region of the global 
minimum. In this paper we are taking the second ap- 
proach and use a grid-search to obtain an approximate 
initial guess. If we examine equations 2, 3, and 1 we 
note that there are three parameters for the binary, r 
(brightness ratio), a, and /3 (separation along the u— 



and w-axes respectively). In addition, there are three 
parameters for each baseline, w, a, and 0. For N base- 
lines there are thus 3x (N + 1) parameters, and it is not 
feasible to do a grid search over that many parameters. 
Fortunately, the coherent integration process results in 
an average phase which is close to zero. If we therefore 
use a binary phase model which has approximately zero 
mean phase (we achieve this by adjusting the vacuum 
path) then we can in many cases reduce the grid search 
to cover just the parameters r, a, and /3. If we make a 
reasonable guess for r from the amplitude of the phase 
variations then the grid search is reduced to just two 
parameters. 

4.1. ^^Tauri 

In the analysis of 9^ Tauri we begin with a grid search 
as outlined above. Figure 4 shows the goodness of fit as 
a function of the two separation parameters. We chose 
a brightness ratio of r = 0.35 based on examining the 
amplitude of the phase oscillations in the raw data. The 
grayscale of the figure indicates the goodness of fit ac- 
cording to a x^ metric, with a brighter color indicating 
a better fit and a darker color indicating a worse fit. 
The best fit, at a « 20 mas, and /3 « 20 mas is framed 
by two white diamonds. Notice that there is also a good 
fit (but not quite as good) at approximately (—a, — /3). 
Finally, notice the complexity of the fitness landscape. 
Without a good initial guess or a optimization method 
which is able to ignore local minima it would be unlikely 
that the fit converges on the global minimum. 

Using this initial guess we then fit the full function 
(Equation 2 with equations 1, and 3 inserted) to the 12 
baseline measurements simultaneously (i.e. using the 
same binary star parameters for every baseline). How- 
ever, we make one modification in that we allow r to 
vary linearly with wavelength. We then obtain the fit 
shown in figure 5. In this figure the solid curve shows 
the baseline phase (including error bars which are often 
too small to distinguish), and the dotted curve shows 
the best-fit model. 

In order to estimate the uncertainty on the fitted 
parameters, r, |j, a, and /3, we perform a bootstrap 
Monte Carlo analysis. Bootstrap error analysis is a sim- 
ple way of estimating the sampling uncertainty from a 
data set with the same distribution as the measured 
data set. It proceeds by repeatedly selecting A^ (with 
replacement) from the N independent segment of the 
data set and performing the full analysis to obtain pa- 
rameters, each time obtaining a slightly different set of 
parameters. If the data set is a good representation of 
the distribution from which it is taken (this is the usual 
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Figure 4. Result of grid scarcli over a, (i on Q Tauri. 

assumption) then the distribution of parameters repre- 
sents the uncertainty of such a sample. We took 100 
bootstrap samples and arrived at separation parame- 
ters of a = 17.134 ± 0.033 mas, /3 = 17.989 ± 0.022 mas, 
which yields a precision on the separation of approx- 
imately 0.16%. The brightness ratio as a function of 
wavelength is plotted in Figure 6. In that figure the 
solid curve is the brightness ratio, whereas the dotted 
curves represent one standard deviation. Also plotted 
are measurements of the brightness ratio from Arm- 
strong et al. (2006)."^ That paper used a much larger 
database than the present work, spanning both Mark 
III and NPOI data, yet arrived at uncertainties in the 
brightness ratio which are similar to or greater than 
ours. Table 2 lists the brightness ratio at several wave- 
lengths. The discrepancy between this work and prior 
work will elaborated on in the discussion section. For 
now we just note that the reduced y^ of the fit is ap- 
proximately eight, corresponding to a mean difference 
between model and data of 2.5 standard deviations. 

4.2. KUrsa Majoris 

The second example, k Ursa Majoris is more challeng- 
ing. It consists of five observations on the same set of 
baselines. We approach this data set in the same way 
as the previous data set, and the final fit to the data is 
presented in Figure 7. The fit corresponds to a bright- 
ness ratio of AM « 0.5, and separation a = —115 mas, 
P = 85 mas. Because of discrepancies between the 
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Figure 5. Visibility phase (solid) and model (dotted) for 6^ Tauri. 
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Figure 6. Brightness ratio as a function of wavelength (solid 
curve) and one standard deviation (dotted curves). Also plotted 
are individual data points taken from."^ 

model and the data we do not specify uncertainties but 
instead refer to the discussion section. 



5. DISCUSSION 

An important consideration when doing high-accuracy 
work is knowledge of the wavelengths to a high degree 



0.50 




0.55 


1.03 ±0.03 


0.65 


1.08 ±0.02 


0.70 


1.10 ±0.01 


0.75 


1.13 ±0.01 


0.80 


1.16 ±0.01 



A (fim) This work Armstrong et al. (2006) 

1.05 ±0.05 
1.12 ±0.03 

1.16 ±0.03 

1.14 ±0.05 



Table 2. Brightness ratio as a function of wavelength for this 
work compared with Armstrong et al. (2006) 



of accuracy. In the case of 9^ Tauri wavelength knowl- 
edge is less critical for estimating the brightness ratio, 
but still important for the separation. In the case of a 
large-separation binary system such as KUrsa Major it 
is much more important. For both the data sets ana- 
lyzed here only nominal wavelength and bandpass infor- 
mation was available, and that likely is a contributing 
factor to the imperfect fit x^ = 8 in the case of 0^ Tauri 
and is probably the dominant effect of the goodness of 
the fit in the case of k Ursa Major. In Figure 7, pan- 
els (b) and (d) contain the same baseline, but recorded 
on two different spectrographs. However the two ob- 
servations are not identical, which is especially evident 
at the blue end of the spectrum. The likely cause of 
this is that the two spectrographs have slightly different 
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Figure 7. Visibility phase (solid) and model (dotted) for k Ursa Majoris. 



wavelength scales. We are not able to correct this effect 
because no measured wavelength information was avail- 
able for that particular date. However on later dates the 
bandpass is measured to a high degree of accuracy. 

Another important thing to notice, if we look at 
Figure 6, is the discrepancy between the present re- 
sults and those of Armstrong ct al. (2006)."^ First, in 
the older data we observe what appears to be a cur- 
vature as a function of wavelength. We do not see this 
in our data because we modeled the brightness ratio to 
be linear with wavelength. It would be interesting to 
add a quadratic term and then repeat the comparison. 



Another thing to note is that Armstrong et al. (2006) 
obtain slightly larger brightness ratios than us. It is 
important to realize in this context that the amplitude 
of the phase oscillation is modulated not only by the 
brightness ratio, but also by the visibility amplitude of 
the component stars. The component stars are unre- 
solved such that visibilities are close to one. We did not 
take that effect into account in this initial work. 

6. CONCLUSION 

We have demonstrated that coherently integrated vis- 
ibility phases can be used to obtain high-quality mea- 



surements of the fundamental parameters of binary 
stars. The measurements are accurate to the point 
where careful calibration of wavelengths and band- 
passes is necessary. This calibration was not available 
for the data sets considered in this paper, but is rou- 
tinely performed on more recent NPOI data sets. 
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